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ABSTRACT 



A FORTRAN program has been developed to allow for the 
use of probabilistic design methods in the numerical optimi- 
zation process. The program was written as a set of sub- 
routines for COPES (Control Program For Engineering Synthesis) 
COPES maximizes or minimizes a numerically defined objective 
function subject to a set of inequality constraints using 
the optimization program CONMIN (A Fortran Program for Con- 
strained Function Minimization) . The program developed here 
allows for the use of both the normal and lognormal distribu- 
tion models. Design examples are presented to demonstrate 
the program capabilities. User instructions are provided for 
inclusion in the COPES user's manual. 
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I. INTRODUCTION 



In the traditional deterministic approach in engineering 
design the variables are treated as single valued numbers. 
These variables, whether dimensions, material properties, 
loads, etc., are actually statistical in nature. By the use 
of safety factors designers protect against these variations 
by usually over designing. A more logical approach is to 
take into account the statistical data known about each 
variable and design for a certain reliability. By consider- 
ing the statistical nature of each variable we should be 
able to better predict reliability and performance. 

Numerical optimization has proven to be a very powerful 
tool in engineering design. Virtually all design problems 
require minimization or maximization of some objective. For 
the design to be acceptable, it must also satisfy a certain 
set of specified requirements called constraints. If these 
constraints are specified as probabilities of failure the 
result would be a probabilistic design. 

The purpose of this research was to test the applica- 
bility of combining probabilistic design concepts with those 
of numerical optimization by: (1) developing a pilot computer 

code to calculate probabilities of failure, (2) incorporate 
this computer code with COPES/CONMIN in order to perform 
probabilistic optimization, and (3) test using numerous 
examples . 
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It is assumed throughout this discussion that the reader 
is familiar with the use of the COPES/CONMIN optimization 
program [Ref. 1], and [Ref. 2]. 
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II. NUMERICAL OPTIMIZATION 



The numerical optimization problem considered here is 
stated as follows: Find the set of n design variables 

contained in the vector X which will 

Minimize F(X) (1) 

Subject to; 

gj(X) < 0 j = l,m (2) 

x^<x. <xV i=l,n (3) 

The components, x^, of X are referred to as design variables 
which are changed to improve the design. The function F(X) 
is called the objective. Inequality constraints, 9 j (X) , are 
the response limits imposed on the design. There are n de- 
sign variables and m inequality constraints. The lower and 

£ u 

upper bounds, x^ and x^ are limits imposed on the design 
variables to insure a practical result. 

For a deterministic design COPES [Ref. 1], determines 
the constraints and objective as follows. Consider the 
design of a single bar undergoing uniaxial tension. If we 
wish to minimize the weight, the objective, 

F(X) = pAL (4) 
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where 



P 



the specific weight; 



A = the area of the bar; 



L 



the length of the bar. 



The design variable is the area, A. If the stress in the 
bar is S and the stress limit imposed is Sy we desire that 
S be less than or equal to Sy so the constraint becomes 



In this investigation a computer subroutine was developed 
to provide constraints based on probabilities of failure. 

If the allowed probability of failure is PF and the 
actual probability of failure was calculated to be Pf, then 
the constraint would be 



S/Sy - 1 < 0 



(5) 



Pf/PF - 1 < 0 



( 6 ) 
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III. PROBABILISTIC DESIGN 



A. RANDOM VARIABLES 

If one has a large population of n elements for which 
some parameter x^,i=l,n is determined, the mean value of 
the population is 



X 



1 ^ 
= - I 



i=l 



(X.) 



The variance is 



V 

X 



1 ^ 
- I 



i=l 






The standard deviation is 



( 7 ) 



( 8 ) 



a 



X 




( 9 ) 



The standard deviation is a measure of the amount of varia- 
bility of the population data. An additional measure of 
the variability is the dimensionless coefficient of variation. 



C 

X 




( 10 ) 



In engineering design the coefficient of variation is fre- 
quently used. 
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B. FUNCTIONS OF SEVERAL VARIABLES 



If y is a function of n variables, the standard deviation 
of the function Y can be approximated as follows ; 

If Y = f(x^; i = l,n) 

Assuming that all x^ are independent random variables 



then according to 


[Ref. 3, p. 59] 


'^Y 


1=1 1 1 



C. STATISTICAL MODELS 

1 . Normal Distribution 

The normal distribution model is one of the most 



widely used. The 


normal density function is given by 


= 


1 «(x-y ) 

^ exp[ . ^ ] (12) 

2a^ 



where : 



fx(x) = 


the probability density function of 
the random variable X; 


'^x 


standard deviation of the random 
variable X; 


^x 


mean of the random variable X. 
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The magnitude of the standard deviation determines the dis- 
persion of the distribution. A large standard deviation will 
result in a wide bell curve and a small standard deviation 
will result in a narrow one. The curve is symmetric about 
the mean. Approximately 68.26% of all samples will lie 
within plus or minus one standard deviation of the mean value. 
The standardized form of the normal density function [Ref. 4; 
p. 194], can be expressed as 



^ exp[-^] 



(13) 



where 



y = 




(14) 



$ (y) = the standardized form of the single 

^ variate normal density function. 

The single variate cumulative normal distribution function 
is given by 



X 

Fj^(x) = / f^(x)dx 



(15) 



where : 



F (x) = the probability that the random variable 

^ X is equal to or less than the specific 

value X. 
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Or, using the expression, 



y 



X - u 

a 



Equation (15) can be rewritten as, [Ref. 4: p. 195] 



Fjj(x) 



1 

/2tt 



(x-Ux)/Oj; 

/ 



2 

exp [^^^] dy 



(16) 



or 



Fx(x) 



$ 

Y 



X - y 



X, 



o. 



(17) 



2 . Lognormal Distribution 

If ln(X) has a normal distribution then X is said to 
have a lognormal distribution. The lognormal density func- 
tion, [Ref. 4: p. 196], is given by 



fx(x) 



/27ro 






-] ; X > 0 



(18) 



X 



2a 



X 



where : 

X = mean of ln(X) ; 
a = standard deviation of ln(X) . 

X 
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The cuinulative lognormal distribution function, [Ref. 4: 
p. 196], is given by 



P^(K) 



ln(x) -X 

1 _ 2 
/ exp[-^]dy; X > 0 

/2tt -00 ^ 



(19) 



or 






X 



The mean, y , and the variance V , of X are 

X X 



( 20 ) 



y^ = exp[X + -^] 



( 21 ) 



V = exp [2X + a^] [exp(aj) - 1] 

•2V XX 



( 22 ) 



Solving for X and a 



X 






(23) 



X 



'-ln(y^) + In [V^ + y^] 



(24) 



D. PROBABILITY OF FAILURE 

Consider the constant area bar in Figure 3.1 under load 
P. If the stress is the bar is S and the stress limit in 
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Figure 3.1. BAR UNDER UNIAXIAL LOAD 

tension for the material is Sy, then the probability of 
failure of the bar is the probability that the stress S will 
be greater than the limit Sy. The method for determining 
the probability of failure in the bar for both the normal 
and the lognormal formats is presented below. 

1 . Normal Distribution 

If both S and Sy are random variables which follow a 
normal distribution then the probability of failure is the 
probability the S is greater than Sy. The failure function, 
[Ref. 5: p. 8-24] 



Z = (Sy - S) 



(25) 



20 



is also normally distributed. The probability of failure 
is the probability that the failure function is less than 
or equal to 0.0 or 



Pf = P(Z< 0) 



(26) 



In the standard form 



Pf 



* f — - 



(27) 



Where, by the algebra of functions, [Ref. 3] 



y 



z 



Sy 



- 



(28) 



and 



a 



z 




(29) 



The safety index, [Ref. 5; p. 8-24], is defined as 



B 




(30) 



The probability of failure is then 



Pf = $g[-B] 



(31) 
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or the area under the unit normal density function from 3 



to 00 . 



pf = f -±z exp[:i|-]d3 (32) 

3 /27t 

2 . Lognormal Distribution 

If the same two variables S and Sy follow a logri- 
thetic normal distribution then the probability of failure 
is expressed by, (Ref. 5: p. 8-27] 

Pf = P(Sy<S) = P(ln(Sy/S) < 0) (33) 



Letting 



Z = ln(Sy/S) 



(34) 



then 



Pf = P(Z < 0) 



(35) 



or 



Pf = <&g[-3] 

The safety index, [Ref. 5: p. 8-29] is 



( 36 ) 
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3 



ln(Sy/S) 



( 37 ) 



Jlnia + cl^) (IH-C^) ] 
where the Sy and S are median values and 

( 38 ) 

( 39 ) 

where 

( 40 ) 

( 41 ) 

3 will be normally distributed and the probability of 
failure is 




Sy = 



Sy 



1 + c 



Sy 



S = 



'1 + c! 



Pf = 



[ 



-ln(Sy/S) 

/ln[(l+C^^) (1+Cg) ] 



] 



( 42 ) 



The probability that S is greater than Sy is then determined 
in the same manner as the normally distributed case, by 
integrating the unit normal density function from 3 to “ . 
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IV. PROBABILISTIC OPTIMIZATION 



The objective of probabilistic optimization is to mini- 
mize or maximize the mean value of a given function subject 
to a set of constraints based on allowed probabilities of 
failure . 

A. CONSTRAINTS 

The computer subroutine developed in this investigation, 
C0PE19, calculates the constraints as described in Chapter 
II. Given a constrained variable in terms of a desired 
probability of failure, one needs only to calculate the 
actual probability of failure to determine the constraint 
value . 

Consider the cantilevered beam in Figure 4.1. A proba- 
bility of failure, PF, is desired for stress in bending. 

The mean values and coefficients of variation for the yield 
strength, Sy, the length L, the load P, and the dimensions B 
and H are given. The stress in bending 



SB 



6PL 



(43) 



The standard deviation of the bending stress 



a 



SB 




(44) 
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Figure 4.1. CANTILEVERED BEAM 
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l<hi-t>| 



The safety index 



3 



gy - SB 



(45) 




The calculated probability of failure 



Pf 




(46) 



The constraint will be 



G 




(47) 



In the computer program developed, C0PE19, the required 
partial derivatives of the constrained variables with 
respect to the variables in the expression are calculated 
by finite difference methods. 

B. OPTIMIZATION 

Consider the three bar truss in Figure 4.2. The geometry 
is specified and there are two independent load conditions. 
The design task is to determine the areas required for each 
bar that will yield a minimum mean value of the structure 
weight while keeping the stress in any member less than the 
yield stress. 
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THREE BAR TRUSS 




P2 \ PI 



Figure 4.2. THREE BAR TRUSS 
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Minimize; W 



pH(/2Al + A2 + /2A3) 



Subject to; Sy compression £ tension 

where SIG^^ is the stress in member i under load condition j. 
Given; Material p = .1 Ib./cu. in. 

Geometry H = 10.0 in. 

Loads PI = P2 = 20,000 lb. 

Stress limits; Sy compression = -15,000 psi 

Sy tension = 20,000 psi 

Maintain symmetry Al = A3 

Beginning with the design Al = A2 = A3 = 1.0 sq. in. and 
optimizing using COPES/CONMIN (deterministic) the following 
results are obtained; 

W = 2.632 lb. 

Al = A3 = .7796 sq. in. 

A2 = .4275 sq. in. 

The critical constraints are SIGH and SIG32 , which are 
equal. The design space is shown in Figure 4.3. As no 
safety factor was used, the above deterministic design 
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f 




■ CONSTANT OBJECTIVE (LBS) 
CONSTANT STRESS SIGH (PS I) 

Figure 4.3. THREE BAR TRUSS DESIGN SPACE (DETERMINISTIC) 
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corresponds to a probability of failure of 0.5 or a 50% 
chance that SIGH and SIG32 will be greater than the stress 
limit Sy. Utilizing the subroutine COPE19 to provide proba- 
bilistic based constraints and assuming that the coefficients 
of variation of all the variables are 0.1, the following 
results are obtained when designing for a maximum probability 
of failure of 0.5 in each member; assuming a normal distribution, 

W = 2.639 lb. 

Al = A3 = .7883 sq. in. 

A2 = .4094 sq. in. 

The design space is shown in Figure 4.4. The lines of con- 
stant probability of failure for SIGH are shown. As one 
would expect the line of constant 0.5 probability of failure 
for SIGH in Figure 4.4 is the same as the constant 20,000 
psi line for SIGH in Figure 4.3. The reverse is also true. 

If the three bar truss is optimized for a maximum probability 
of failure of 0.4 the resulting safety factor for SIGH is 
1.03 or 19,361 psi and the minimum weight is 2.726 lb. Using 
19,361 psi as the limit of stress for a deterministic optimi- 
zation yields approxiamtely the same minimum weight. In both 
cases the feasible design space is the same with the line 
of constant 0.4 probability of failure and constant 19,361 
psi lying on top of each other. 
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CONSTANT OBJECTIVE (LBS) 

CONSTANT PROBABILITY OF 

FAILURE SIGH 

Figure 4.4. THREE BAR TRUSS DESIGN SPACE (PROBABILISTIC) 
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Table I shows the results of minimizing the weight of 
the three bar truss for various probabilities of failure. 
Table II provides the deterministic optimum designs for the 
equivalent factor of safety. The difference in designs is 
within the numerical accuracy of the optimization. 



TABLE I 

THREE BAR TRUSS PROBABILISTIC DESIGNS 



PROBABILITY 


WEIGHT 


Al = A3 


A2 


SIGH 


SAFETY 


OF FAILURE 


LBS. 


SQ. IN. 


SQ. IN. 


PSI 


FACTOR 


.5 


2.6390 


.78827 


.40938 


20000 


1.00 


.4 


2.7260 


. 81606 


.41785 


19361 


1.03 


. 3 


2. 8238 


. 84478 


.43438 


18691 


1.07 


.2 


2.9430 


.86887 


.48544 


17940 


1.11 


.1 


3.1208 


.92482 


.50500 


16914 


1.18 


.01 


3.6099 


1.0490 


.64280 


14640 


1.37 


.001 


4.0459 


1.1777 


.71493 


13060 


1.53 


.0001 


4.4770 


1. 3212 


.73994 


11792 


1.70 





TABLE II 

THREE BAR TRUSS DETERMINISTIC 


DESIGNS 


SAFETY 


WEIGHT 


Al = A3 


A2 


SIGH 


FACTOR 


LBS. 


SQ. IN. 


SQ. IN 


PSI 


1.00 


2.6326 


. 77962 


.42752 


20051 


1.03 


2.7212 


. 80145 


.45437 


19402 


1.07 


2.8200 


.82773 


.47889 


18726 


1.11 


2.9322 


. 87045 


.47014 


18002 


1.18 


3.1180 


.91845 


.52000 


16933 


1.37 


3.5936 


1.0807 


.53693 


14688 


1.53 


4.0258 


1.2033 


.62227 


13110 


1.70 


4.4800 


1.3379 


.69573 


11781 


It 


should be 


noted that the correspondence of the proba- 


bility 


of failure 


and the safety 


factor 


shown here is unique 
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to this problem, in which all C = 0.1. In general nonlinear 
optimization, this one to one relationship cannot be assured. 

C. EFFECTS OF COEFFICIENTS OF VARIATION ON THE DESIGN SPACE 

As the coefficients of variation of the variables used in 
the design subroutine become smaller, the bands in the design 
space for the constrained values between Pf = 1.0 and Pf = 0.0 
becomes narrower. Figure 4.5 shows the design space for the 
three bar truss in the vicinity of the optimum for the case 
where all variables have coefficients of variation, C = 0.1. 
Figure 4.6 shows the same space for the case where all coeffi- 
cients of variation, C = 0.001. In both figures the lines of 
constant objective and probabilities of failure for the 
critical constrained variable, SIGH, are shown. The line 
of constant Pf = 0.5 is the same for both but the band between 
Pf = 0.9 and Pf = 0.1 is much narrower for the case where 
C = 0.001. 

This phenomenon can cause numerical difficulties in the 
optimization process. Premature termination of optimization 
can occur as the optimizer may be unable to move down the 
resulting 'narrow valley'. The absolute and relative termina- 
tion criteria of CONMIN may have been satisfied but the true 
optimiim might not have been reached. 

If the initial design is a considerable distance from 
the narrow band of changing probabilities of failure, CONMIN 
may obtain zero gradients for the constraints which would 
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R2 (SQ. IN 




CONSTANT OBJECTIVE (LBS) 



CONSTANT PROBABILITY OF FAILURE 
SIGH 



Figure 4.5. THREE BAR TRUSS DESIGN SPACE (C = 0.1) 
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R2 (SQ. IN 




— CONSTANT OBJECTIVE (LBS) 

CONSTANT PROBABILITY OF FAILURE 

SIGH 

Figure 4.6. THREE BAR TRUSS DESIGN SPACE (C = 0.001) 



35 



also cause termination of the optimization process . A 
routine to attempt to prevent this was included in C0PE19 . 

The cumulative distribution curve for the normal format is 
shown in Figure 4.7. The numerical integration routine used 
to calculate the area under the unit normal distribution 
curve yields a probability of failure, Pf of 1.0000000 for 
a S of -5.2999821 and a Pf of 0.0000000 for a 6 of 5.5000086. 
It was modified to "widen" the zone between Pf = 1.0 and 
Pf = 0 . A polynomial (a line equation) was placed at each 
end of the cumulative distribution routine to yield a Pf 
of 1.0 at 6 equals -25.0, decreasing to a Pf of 0.9999866 
at -4.2 and a Pf of 0.00001335144 at g equals 4.2 decreasing 
to a Pf of 0.0 at 25.0. This is shown in Figure 4.8. This 
provides the same optimization results within the accuracy 
of the program, but avoids the problem of zero gradients of 
the constrained variables. 

An additional approximation was used to determine the 
probability of failure when the coefficients of variation 
concerned are zero or nearly so. The probability that X is 
greater than Y is, in the standard format 





(48) 



where 



(49) 



z 



y 



X 
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Figure 4.7. CUMULATIVE NORMAL DISTRIBUTION CURVE 
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Figure 4.8. MODIFIED CUMULATIVE NORMAL DISTRIBUTION CURVE 
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and 



a 

2 




+ a 



(50) 



Mathematically , if a — 0 , then one of three situations can 
occur. 



1. 


X = 


Y, 


6 = 


0, then Pf = 0.5; 


2. 


X > 


Y, 


B ^ 


-<», then Pf = 1.0; 


3. 


X < 


Y, 


6 ^ 


<», then Pf = 0.0. 



This is shown graphically in Figure 4.9. In order to provide 
CONMIN with a smooth function with which to work, a polynomial 
was used to connect the Pf = 1.0 and Pf = 0.0 in Figure 4.9. 
This is shown in Figure 4.10. The interpolating polynomial 
is used when Q, the modified safety index, is within one 
unit of 0. Where S is defined to be 

6 = - (51) 

A similar approximation was developed for the lognormal 
model. 

Experience has shown that these modifications dramatically 
improve the numerical stability of the optimization process 
without significantly affecting the mathematical result. 
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Figure 4.9. CUMULATIVE NORMAL DISTRIBUTION CURVE 
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V. COMPUTER PROGRAM 



A. COPES INPUT DATA 

Standard deterministic optimization requires that COPES 
data blocks A—0 and V be used. For probabilistic optimiza- 
tion the same data is entered in these blocks with the 
exception of NCALC and IPROB in data block B. Here an NCALC 
= 7 would be entered to indicate probabilistic optimization 
was to be performed and an IPROB = 1 or 2 to indicate which 
probability model is to be used. Additional data blocks are 
required for probabilistic optimization. Appendix A con- 
tains detailed input instructions for these data blocks as 
well as an example problem including sample input data. 

B. OPTIMIZATION DATA FLOW 

The interaction between COPES, CONMIN , and the analysis 
subroutine will be the same for probabilistic optimization 
as it is in standard deterministic optimization with one 
exception. In the optimization process, each time COPES 
requires the determination of a constraint vector it will 
call COPE 19. C0PE19 will call the analysis subroutine as 
many times as necessary to determine the required partial 
derivatives by finite difference steps. Once these are 
obtained, the constraint vector is determined and provided 
to COPES. A flow diagram for probabilistic optimization 
is provided in Figure 5.1. 
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Figure 5.1. PROBABILISTIC OPTIMIZATION FLOW DIAGRAM 



C. COPE 19 SUBROUTINE 

The main computer subroutine developed in this investi- 
gation is C0PE19. Three additional subroutines, COPE20 
through COPE22, were also developed to perform specific 
tasks required by C0PE19 . The operation and flow of these 
subroutines are discussed below. 

1. Input 

The following is provided to C0PE19: 

a. The COPES control parameter ICALC. 

b. Which probability model to use, either normal 
or lognormal. 

c. The number of probability variables which make 
up the constraint equations (IVAR) . 

d. The global locations in the ANALIZ subroutine 
common block of each IVAR. 

e. Coefficient of variation of each IVAR, assumed 
constant throughout the optimization. 

f. The number of constrained variables. 

g. The global location in the ANALIZ subroutine 
common block of each of the constrained 
variables . 

h. The upper and lower limits imposed on each of 
the constrained variables. 

i. The allowed probability of failure at the upper 
and lower limits of each constrained variable. 

j. The coefficient of variation of the constrained 
variable limits . 
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k. The variable MGRAD which determines how often 

the constrained variable gradients are calculated. 

2 . Determination of Constraint Vector 

When COPES requires a constraint vector for a given 
design it will call C0PE19 with ICALC = 2. C0PE19 will then 
calculate the partial derivatives of each constrained varia- 
ble with respect to the variables which make up the con- 
straint equations. The subroutine COPE20 controls the finite 
difference calculations. In the example in Section IV. A, 
the cantilevered beam, it would calculate: 



9SB 

3L 



3SB 

3P 



3SB 

3H 



3SB 

3B 



6P 

br2 

6L 



12PL 

6PL 

2 2 
B H 



(52) 

(53) 

(54) 

(55) 



C0PE19 will then calculate the standard deviation (Equation 
44) of the constrained variable, in this case SB, as dis- 
cussed in Section III.B and Section IV. A. Now the safety 
index is calculated and subroutine C0PE21 called to deter- 
mine the probability of failure. When required, C0PE21 will 
call subroutine COPE22 to perform numerical integration of 
the normal density function. The constraint value to be 
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stored in the G vector is then determined (Equation 6) . 

The above would be performed for each constrained variable 
and the resulting G vector provided to COPES. 

D. REDUCTION OF FUNCTION EVALUATIONS 

In order to determine the gradients of the constrained 
variables finite difference steps are taken resulting in a 
significant increase in the number of analyses performed in 
the optimization process. In many problems these gradients 
remain essentially the same throughout the optimization. 
Therefore, it is sometimes possible to reduce the frequency 
with which the gradients are calculated and achieve approxi- 
mately the same result. The input variable code LGRAD is 
included as an option in probabilistic optimization to allow 
this choice. LGRAD has the following meanings: 

0: Calculate gradients each time COPE19 is called. 

1; Calculate gradients at the beginning of each 
CONMIN iteration. 

2; Calculate gradients only at the beginning of 
optimization. 

Table III demonstrates the result of various LGRAD selections 
when designing the three bar truss for an allowed probability 
of failure of 0.50, where all coefficients of variations, 

C = 0.1. 

Probabilistic optimization was performed using the three 
LGRAD values on numerous additional test cases. In almost 
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TABLE III 



REDUCTION OF FUNCTION EVALUATIONS: THREE BAR TRUSS DESIGN 





LGRAD 


WEIGHT 


A1 = A3 


A2 


NUMBER OF 






LBS. 


SQ. IN. 


SQ. IN. 


FU^'CT EVALS 


Deterministic 




2.6326 


.77962 


.42752 


35 


Probabilistic 


0 


2.6390 


. 78827 


.40938 


343 


Probabilistic 


1 


2.6389 


. 77386 


.45011 


244 


Probabilistic 


2 


2.6379 


. 78972 


.40427 
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every case the optimization results were essentially the same 
for LGRAD = 0 and 1. The composite driveshaft design in 
Chapter V provides an excellent example showing the function 
evaluation reductions when performing a relatively complex 
design. Table IV demonstrates the results of reducing the 
frequency of constrained variable gradient calculations 
when designing the composite driveshaft for minimum weight, 
using the normal distribution model. For this particular 
problem, calculation of constrained variable gradients only 
at the beginning of optimization was insufficient to obtain 
results near the optimum. 



TABLE IV 

REDUCTION OF FLT^CTION EVALUATIONS: DRIVESHAFT DESIGN 



LGRAD WEIGHT 

LBS . 

0 9.0203 

1 9.0205 

2 27.216 



NUMBER OF FUNCTION 
EVALUATIONS 

4011 

1344 

89 
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VI . NUMERICAL EXAMPLES 



The following examples are presented to demonstrate the 
capabilities of the computer subroutines developed in this 
investigation. In both cases an existing ANALIZ subroutine 
was used after minor modifications. For these two shaft 
designs the subroutine from [Ref. 6] , was used. The ANALIZ 
subroutine was altered to include the required probability 
variables in the global common block and to remove the 
factor of safety calculations which were originally used as 
constraints . 

In each example a similar deterministic design was per- 
formed to give the reader a better feel for the results. 

The deterministic designs are not intended to duplicate the 
probabilistic design results. 

The mean values used in the designs are from [Ref. 6] . 
The coefficients of variation are reasonable assumed values 
used for demonstration purposes. 

A . STEEL DRIVESHAFT 

Design, for minimum weight, a steel driveshaft to trans- 
mit 150.0 horsepower at 300.0 RPM. The loading is presented 
in Figure 6.1. The shaft is designed against failure in 
strain, in torsional buckling, vibration frequency and a 
maximum deflection of 0.05 inches. 

The design variables are the shaft thickness and the 
shaft inside diameter. The thickness is to be maintained 
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F = 50 LB 
T= 31500 IN LB 



Figure 6.1. DRIVESHAFT LOADING 
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between 0.01 and 2.0 inches and the inside diameter between 
0*5 and 5.0 inches. The initial design is an inside 

diameter of 2.0 inches and a thickness of 1.0 inch. The 
remaining data for the analysis is presented in Table V. 



TABLE V 

INPUT DATA FOR STEEL DRIVESHAFT 



HORSEPOWER 

SPEED 

AXIAL LOAD 
LENGTH 

YOUNGS MODULUS 
SHEAR MODULUS 
POISSONS RATIO 
SPECIFIC WEIGHT 
STRAIN LIMIT 
MAX DEFLECTION 
THICKNESS 
INSIDE DIAMETER 
FAILURE DUE TO: 

TORSIONAL BUCKLING 
FREQUENCY 



MEAN 

VALUE 

150. HP 
300. RPM 
50. LB. 
120. IN. 
30,000,000 PSI 
11,538,000 PSI 
.3 

.282 LB/CU. 
.001 IN. /IN 
.05 IN. 



COEFFICIENT 
OF VARIATION 

.08 

.08 

.18 

.015 

.038 

.042 

.026 

IN. .01 

.01 
.05 
.015 
.015 



.21 

.06 



Allowing for maximum probabilities of failure of 0.01 in 
strain, torsional buckling, deflection and vibration fre- 
quency; a minimum weight of 54.15 lbs is achieved for the 
normal model and 52.88 lbs for the lognormal model. The 
initial design summary is presented in Table VI and the 
final designs summaries and COPES optimization results in 
Tables VII through XII. 

Designing for a factor of safety against failure of 2.0 
for strain, torsional buckling, vibration frequency, and 
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TABLE VI 



STEEL DRIVESHAFT INITIAL DESIGN SDMMARY 



STEEL DRIVESHAFT OUTPUT 



NUMBER OF PLYS = 1 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = D.O 



DIMENSIONS 

PLY THICKNESS PERCENT DIAMETER THETA 

INSIDE 0.20000E+01 

1 0. 10000E+01 100.00 0.40000E+01 0.0 



STIFFNESS 
AE = 0.31071E+09 
El = 0.38838E+09 
GJ = 0.27187E+09 



LOADS : 

L C T 

i'"*0.31513E+05 0‘.‘0 



F H 

0.50000E+02 0.1505oe+03 0. 



LOAD CONDITION 1 

PTY PPT P PP'T’ 

1 0.16092E-06 lOOloO oVo 

CRITICAL SPEED 
MAXIMUM DEFLECTION = 
WEIGHT 
VOLUME 



S.F. 

100.00 

0 . 1554 3E+0'4 
0 . 19195E-01 
0 .31893E+03 
0 . 1 1310E+04 



EPLT 

0. 23182E-03 



RPM 

30000E+03 



S.F. 

5.61 
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TABLE 7II 



STEEL DHIVESHAFT COPES STANDARD OUTPOT: NORHAL 



OPTIMIZATION RESULTS 



OBJECTIVE FUNCTION 

GLOBAL LOCATION 3 FUNCTION VALUE 0.54154E+02 



DESIGN VARIABLES 

D. V. GLOBAL 
ID NO. VAR. NO. 
1 1 11 

2 2 2 



DESIGN CONSTRAINTS 

DETERMINISTIC 





GLOBAL 


LOWER 


MEAN 


UPPER 


ID 


VAR. NO. 


BOUND 


VALUE 


SOUND 


1 


51 


-0. 10000S-0 2 


0.94774E-06 


0. 10000E-02 


3 


52 


-0. 10000E-0 2 


0.0 


0. 10000E-02 


5 


53 


-0. 10000 E-0 2 


0.6821 3E-03 


0.10000E-02 


7 


54 


0 . 31 513E+0 5 


0.24340E+06 


0. 1 1000E+16 


8 


56 


0.30000E+03 


0.25063E+04 


0.1 1000E+16 


9 


57 


-0. 11000E+16 


0.72076E-02 


0.50000E-01 



LOWER 

BOUND 

0.10000E-01 
0.5000 OE+00 



VALUE 

0. 99883E-01 
0. 50000E+01 



UPPER 

BOUND 

0.20000E+01 

0.50000E+01 
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ID 

1 

3 

5 

7 

8 

9 

ID 

1 

3 

5 

7 

8 

9 



TABLE VIII 



STEEL DHIVESHAPT COPES PROBABILISTIC OUTPUT: NORMAL 



PROBABILISTIC 



MEAN 

VALUE 

0 .94774S-06 

0.0 

0.68213E-03 

0.24340E+06 

0.25068S+04 

0.72076E-02 



STANDARD 
DEVIATION 
0. 1753S-06 
0.0 

0. 6535E-04 
0. 1364E + 05 
0. 9634E+02 
0.5707E-03 



COEFFICIENT 
OF VARIATION 
0. 1753E-01 
0.0 

0. 9581S-01 
0. 5606E-01 
0. 3843E-01 
0.7918E-01 



PROBABILITY OF FAILURE 
AT LOWER BOUND 



PROBABILITY OF FAILURE 
AT UPPER BOUND 



ALLOWED 
0. 1000E-01 
0. 1C00E-01 
0. 1000E-01 
0. 1000E-01 
0. 1000E-01 
0. 1000E + 01 



CALCULATED 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



ALLOWED 
0. 1000S-01 
0. 1000E-01 
0. 1000E-01 
0. 1000B+01 
0. 1000E + 01 
0. 1000E-01 



CALCULATED 

0.0 

0.0 

0. 1000E-01 
0 . 0 
0.0 
0.0 
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TIBLE IX 



STEEL DRI7ESHAFT DESIGN SOMNARY: NORMAL 



DIMENSIONS 

FLY THICKNESS PERCENT DIAMETER THETA 

INSIDE 0.50000E+01 

1 0. 99883S-01 100,00 0.51993E+01 0.0 



STIFFNESS 
AE = 0,52757E+08 
El = 0.17158E+09 
GJ = 0.12011E+09 



STEEL DRIVESHAFT OUTPUT 



NUMBER OF PLYS 
NUMBER OF MATERIAL TYPES 
NUMBER OF LOAD CONDITIONS 
ECCENTRICITY 



0.0 



LOADS 

L • ^ . 



C T 

r5. 31513E+05 



M 

0 . 0 



H.P. 

0. 15000E+03 



RPM 

0.30000E+ 



0.50 OOOB + 02 



LOAD CONDITION 1 
PLY EPL S.F. 

1 0.94774E-06 100.00 



0.0 



SPT 



S.F. 

100.00 



EPLT 

0. 68213E-03 



S.F. 

1.91 



CRITICAL SPEED 
MAXIMUM DEFLECTION 
WEIGHT 
VOLUME 




0 .54154E+02 
0 . 192 04E+03 



TABLE X 



STEEL DRIVESHAPT COPES STANDARD OOTPUT; LOGNORMAL 



OPTIMIZATION RESULTS 



OBJECTIVE FUNCTION 

GLOBAL LOCATION 3 FUNCTION VALUE 0.52882E+02 



DESIGN VARIABLES 



ID 

1 

2 



D. V. 

NO. 

1 

2 



GLOBAL 
VAR. NO. 
11 
2 



LOWER 

BOUND 

0. 10000E-01 
0.5000 OE+00 



VALUE 

0. 97581E-01 
0. 50000E+01 



UPPER 

BOUND 

0.20000E+01 

0.50000E+01 



DESIGN CONSTRAINTS 

DETERMINISTIC 





GLOBAL 


LOWER 


MEAN 


UPPER 


ID 


VAR. NO. 


BOUN D 


VALUE 


BOUND 


1 


51 


-0. 10000 E-02 


0.97054E-06 


0. 10000E-02 


3 


52 


-0 . 10000 E-02 


0.0 


0. 10000F-02 


5 


53 


-0.10000E-02 


0. 55856S-03 


0. 10000E-02 


7 


54 


0 .31513E+05 


0.22961E+06 


0. 11000E+16 


8 


56 


0.30000E+03 


0. 25056E + 04 


0. 11000E+16 


9 


57 


-0.1 1000 E+16 


0.72144E-02 


0 .50C00E-01 
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TABLE XI 



STEEL DRIVESHAET COPES PROBABILISTIC OUTPUT: LOGNORMAL 



PROBABILISTIC 



MEAN 

VALUE 

0. 97054E-06 

0.0 

0. 69856E-03 
0. 22961E+06 
0. 25056E + 04 
0. 72144E-02 



STANDARD 
DEVIATION 
0. 1796B-06 
0.0 

0.6692E-04 
0. 1283E+05 
0.9630E+02 
0.5713E-03 



COEFFICIENT 
OF VARIATION 
0. 1796E-01 
0.0 

0.9580E-01 
0.561 1S-01 
0.3843E-01 
0.7918E-01 



PROBABILITY OF FAILURE 
AT LOWER BOUND 



ID ALLOWED 

1 0.1000E-01 

3 0.1000E-01 

5 0.1000E-01 

7 0.1000E-01 

8 0.1000E-01 

9 0.1000E+01 



CALCULATED 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



PROBABILITY OF FAILURE 
AT UPPER BOUND 



ALLOWED 
0. 1000E-01 
0. 1000E-01 
0. 1000E-01 
0. 1000E+01 
0. 1000E+01 
0. 1000E-01 



CALCULATED 
0 . 0 
0.0 

0. 9993E-02 
0 . 0 
0 . 0 
0 . 0 
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raBLE XII 



STEEL DRI7ESHAPT DESIGN SOHMARY: LOGNORMAL 



STEEL DRIVESHAPT OUTPUT 



NUMBER OP PLYS = 1 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = 0.0 



DIMENSIONS 

PLY THICKNESS PERCENT DIAMETER THETA 

INSIDE 0.50000E+01 

1 0.97581E-01 100.00 0.51952E+01 0.0 



STIFFNESS 
AE = 0.51518E+08 
El = 0.16740E+09 
GJ = 0.11718E+09 



LOADS : 

L.C. T M 


p 


H. P. 


RPM 


1 0.31513E+05 0.0 


0.50000E+02 


0. 15000E+03 


0.30000E+03 



LOAD CONDITION 1 
PLY SPL S.F. 

1 0. 97054 E-06 100.00 

CRITICAL SPEED 
MAXIMUM DEFLECTION 
WEIGHT 
VOLUME 



EPT 

0.0 



^ 00 \ 00 



= 0.25056E+04 
= 0 .72144E-02 
= 0 .52882E+02 
= 0 . 18752E+03 



TT-OT T S 

0.69856E-03 1 ! 86 
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deflection an optimum of 73.78 lbs is obtained. The COPES 
optimization results and final design summary are contained 
in Tables XIII and XIV. 

B. COMPOSITE DRIVESHAFT 

Design, for minimum weight a 120 in. long, four ply 
graphite epoxy driveshaft to transmit 150.0 horsepower at 
300 RPM. Figure 6.1 shows the shaft loading condition. 

The shaft is designed against failure in transverse, logitu- 
dinal, and shear strain; torsional buckling, vibration fre- 
quency and a maximum deflection of 0.05 inches. Table XV 
contains the initial design. 

The design variables are the thickness of each ply, the 
orientation of each ply and the inside diameter. Addition- 
ally, the thickness of ply 2 and ply 3 must remain equal to 
each other and the orientation of ply 3 equals the negative 
of ply 2. The following additional constraints are placed 
on the design variables; The inside diameter must remain 
between 0.5 inches and 5.0 inches and the thickness of any 
given ply between 0.01 and 0.5 inches. Ply 2 is allowed to 
vary between 0.0 and 90.0 degrees. The orientation of plys 
1 and 4 remain constant. The remaining input data is pre- 
sented in Table XVI. 

Allowing for maximum probabilities of failure of 0.01 in 
strain, torsional buckling, deflection, and vibration fre- 
quency; an optimiom weight of 9.02 lbs is achieved using the 
normal distribution model and 10.28 lbs for the lognormal 
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TABLE XIII 



STEEL DHI7ESHAFT COPES OOTPOT : DETERMINISTIC 



OPTIMIZATION RESULTS 



OBJECTIVE FUNCTION 

GLOBAL LOCATION 3 FUNCTION VALUE 0.73783E+02 



DESIGN VARIABLES 

D. V. GLOBAL 
ID NO. VAR. NO. 
1 1 11 



LOWER 

BOUND VALUE 

0. 1 OOOOE-0 1 0. 1 3515E+00 

0.50000E+00 0.50000E+01 



UPPER 

BOUND 

0.20000E+01 

0.50000E+01 



DESIGN CONSTRAINTS 

DETERMI NISTIC 





GLOBAL 


LOWER 


MEAN 


UPPER 


ID 


VAR. NO. 


BOUN D 


VALUE 


BOUND 


1 


51 


-0.50000E-03 


0.59561E-06 


0.50000E-03 


3 


52 


-0.50000E-03 


0.0 


0. 50000E-03 


5 


53 


-0.50000 E-03 


0.50035E-03 


0 .50000E-03 


7 


54 


0 .63025 E + 05 


0.51838E+06 


0. 11000E+16 


8 


56 


0.60000 E+03 


0. 25247E+04 


0. 11000E+16 


9 


57 


-0.11000 E+16 


0.7 1044E-02 


0. 20000E-01 
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T&BLE XIV 



STEEL DRI7ESHAFT DESIGN SOMMARY; DETERHINISTIC 



STEEL DRIVESHAFT OUTPUT 



NUMBER OF PLYS = 1 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = O.d 



DIMENSIONS 

PLY THICKNESS PERCENT DIAMETER THETA 
INSIDE 0.53000E+01 

1 0. 13515E + 00 100,00 0.52703E+01 0.0 



STIFFNESS 
AE = 0.71879E+08 
El = 0. 237 0SE+09 
GJ = 0.16597E+09 



LOADS : 

L.C. T M 

1 0.31513E+ 05 0.0 



p 

0.50000E+02 



H. P. 

0. 15000E+03 



RPM 

0.30000E+03 



LOAD CONDITION 1 

PTY FPT 9 F 

1 0.69561E-06 100.00 

CRITICAL SPEED 
MAXIMUM DEFLECTION 
WEIGHT 
VOLUME 



EPT S.F. 

0.0 100.00 

= 0 .25247E+0'-; 

= 0 .71044E-02 
= 0 .73783E+02 
= 0 . 26164E+03 



EPLT S.F. 

0.50035E-03 2.60 
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TABLE XV 



COMPOSITE DRIVESHAFT INITIAL DESIGN 



THICKNESS 

IN. 



PLY 1 .25 
PLY 2 .25 
PLY 3 .25 
PLY 4 .25 

INSIDE DIAMETER 2.0 IN. 



ORIENTATION 

DEGREES 

0.0 

20.0 

- 20.0 

90.0 



TABLE XVI 

INPUT DATA FOR COMPOSITE DRIVESHAFT 



MEAN COEFFICIENT 

VALUE OF VARIATION 



HORSEPOWER 

SPEED 

LENGTH 

AXIAL LOAD 

LONGITUDINAL MODULUS 

TRANSVERSE MODULUS 

SHEAR MODULUS 

MAJOR POSSIONS RATIO 

STRAIN LIMITS; 

LONGITUDINAL COMPRESSIVE 
LONGITUDINAL TENSILE 
TRANSVERSE COMPRESSIVE 
TRANSVERSE TENSILE 
MAX SHEAR 
SPECIFIC WEIGHT 
MAX DEFLECTION 
PLY ORIENTATION 
PLY THICKNESS 
INSIDE DIAMETER 
FAILURE DUE TO: 

TORSIONAL BUCKLING 
FREQUENCY 



150 . 


HP 


.08 


300 . 


RPM 


.08 


120. 


IN. 


.015 


50. 


LB. 


.18 


21,000,000 


PSI 


. 14 


1,700,000 


PSI 


.14 


650,000 


PSI 


. 14 


.21 




.08 


-.00857 


IN. /IN. 


.15 


.00857 


IN. /IN. 


.15 


-.0176 


IN. /IN. 


.15 


.00471 


IN. /IN. 


.15 


.0184 


IN. /IN. 


.15 


.056 


LB/CU. IN. 


.05 


.05 


IN. 


.05 

.05 

.015 

.015 

.20 

.06 
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model. The ANALIZ subroutine design summary for the initial 
design is presented in Table XVII. Tables XVIII through 

XXIII contain the design summaries and final COPES optimiza- 
tion results. 

Designing for factors of safety against failure of 2.0 
for strain, torsional buckling, deflection, and vibration 
frequency an optimum weight of 10.23 lbs is obtained. Tables 

XXIV and XXV present the deterministic design results. 
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TABLE XVII 



COHPOSITE DRIVESHAFT INITIAL DESIGN SUMMARY 



COMPOSITE DRIVESHAFT OUTPUT 



NUMBER OF PLYS = '4 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = 0.0 



DIMENSIONS 
PLY THICKNESS 
INSID 

1 0.25000E+00 

2 0. 25000E + 00 

3 0.25000E+00 

4 0. 25000E + 00 



PERCENT DIAMETER 
E 0.20000E+01 

25.00 0.25000E+01 

25.00 0.30000E+01 

25.00 0.35000E+01 

25.00 0.40000E+01 



THETA 

0.0 

0. 20000E+02 
-0.20000E+02 
C.90000E+02 



STIFFNESS 
AE = 0. 12143E+09 
El = 0.12426E+09 
GJ = 0.37225E+08 



LOADS: 



L.C 


• 


T 


M 


F 


H. 


P. RPM 


1 


0 


.31513E+ 05 


0. 0 


0.50000E+02 


0.15000S+03 0.30000E+03 


LOAD 


CONDITION 


1 










PLY 




EPL 


S. F. 


EPT 


S. F. 


EPLT 


S.F. 


1 


0. 


41 176E-06 


100. 00 


0.0 


100. 00 


0.10582E-02 


17.39 


2 


0. 


40847E-03 


20. 98 


-0.40 806E-03 


43.13 


0. 97246E-03 


18.92 


3 - 


0. 


47576 S-03 


18. 01 


0.47617E-03 


9.89 


0. 11351E-02 


16.21 


4 


0. 


53159E-09 


100. 00 


0.41 123E-06 


100.00 


-0. 16931E-02 


10.87 






CRITICAL SPEED 


= 0 . 19725E+04 










MAXIMUM DEFLECTI 


ON = 0.1 1746E- 


01 










WEIGHT 




= 0 .63334E+02 










VOLUME 




= 0 .11310E+04 
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T&BL3 X7III 



COMPOSITE DEIVESHAPT COPES STANDARD OOTPOT: NORMAL 

OPTIMIZATION RESULTS 
OBJECTIVE FUNCTION 

GLOBAL LOCATION 3 FUNCTION VALUE 0,90203E+01 



DESIGN VARIABLES 





D. V. 


GLOBAL 


LO WER 




UPPER 


ID 


NO . 


VAR. NO. 


BOUND 


VALUE 


BOUND 


1 


1 


11 


0. 10000E-01 


0. 10000E-01 


0.50000E+00 


2 


2 


12 


0.10000S-01 


0. 1C212E+00 


0.50000E+00 


3 


2 


13 


0. 10000E-01 


0. 10212E+00 


0.50000E+00 


4 


3 


14 


0.10000E-01 


0. 10000E-01 


0.50000E+00 


5 


4 


32 


0.0 


0. 34968E+02 


0.90000E+02 


6 


4 


33 


0.0 


-0. 34968E+02 


-0.9000CE+02 


7 


5 


2 


0.50000E+00 


0. 1631 1E + 01 


0.50000E+01 



DESIGN 


CONSTRAINTS 












DETERMI NISTIC 








GLOBAL 


LONER 


MEAN 


UPPER 


ID 


VAR. NO, 


BOUND 


VALUE 


BOUND 


1 


51 


-0 .85700 E-02 


0. 35785S-05 


0. 85700E-02 


3 


52 


-0, 17600E-01 


0.0 


0.47 100E-02 


5 


53 


-0 .18400 E-01 


0.47679S-02 


0. 18400E-01 


7 


54 


-0 .85700S-02 


0.25105E-02 


0. 85700S-02 


9 


55 


-0.17600 E-01 


-0, 25069E-02 


0. 47100E-02 


1 1 


56 


-0.18400S-01 


0. 18288E-02 


0. 18400E-01 


13 


57 


-0.85700 E-02 


-0. 27745E-02 


0.85700E-02 


15 


58 


-0 . 17600E-01 


0. 2778 IE-02 


0.47100E-02 


17 


59 


-0 .18400 E-01 


0. 20319E-02 


0. 18400E-01 


19 


60 


-0.85700B-02 


0. 18741E-08 


0. 85700E-02 


21 


6 1 


-0 . 17600 S-01 


0. 35766E-05 


0. 47100E-02 


23 


62 


-0.18400 E-01 


-0, 59688E-02 


0, 18400E-01 


25 


63 


0 .31 51 3E+05 


0.55160E+05 


0. 11000E+16 


26 


65 


0.30000E+03 


0. 1 1715E+04 


0. 11000E+16 


27 


66 


-0.11000 E+ 16 


0. 3481 5E-01 


0.50000E-01 
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TABLE XIX 



COMPOSITE DRIVESHAFT COPES PROBABILISTIC ODTPlJT; NORMAL 



PROBABI LI STIC 





MEAN 


STANDARD 


COEFFICIENT 


ID 


VALUE 


DEVIATION 


OF VARIATION 


1 


0.35785E-05 


0 .4929E-06 


0.4929E-01 


3 


0.0 


0.0 


0.0 


5 


0.47679E-02 


0.7317E-03 


0. 1535E+00 


7 


0.251 05E-02 


0. 3838E-03 


0. 1529E+00 


9 


-0.25069E-02 


0.3834S-03 


0. 1530E+C0 


1 1 


0.18288E-02 


0. 4363E-03 


0. 2386E+00 


13 


-0. 277U5E-02 


0.4221E-03 


0. 152 1E+00 


15 


0.27781E-02 


0.4225E-03 


0. 1521E+00 


17 


0.20319E-02 


0.4938E-03 


0. 2430S+00 


19 


0. 18741E-08 


0.4687E-03 


0. 1000E+02 


21 


0.35766E-05 


0.4687E-03 


0. 1000E+02 


23 


-0.59688E-02 


0.9181E-03 


0. 1538E+00 


25 


0.55160E+05 


0. 7849E+04 


0. 1423E+00 


26 


0 . 1 1715E+04 


0. 9463E+02 


0. 8078E-01 


27 


0.34815E-01 


0.6032E-02 


0. 1733E+00 




PROBABILITY OF FAILURE 


PROBABILITY OF FAILURE 




AT LOWER BOUND 


AT UPPER BOUND 


ID 


ALLOWED CALCULATED 


ALLOWED 


CALCULATED 


1 


O.lOOOE-01 0.0 


0. 1000E-01 


0.0 


3 


0.1000E-01 0.0 


0. 1000E-01 


0. 0 


5 


0.1000E-01 0.0 


0. 1000E-01 


0. 9537E-06 


7 


0.1000E-01 0.0 


0. 1000E-01 


0.3338E-05 


9 


0.1000E-01 0.0 


0. 1000E-01 


0. 0 


1 1 


0.1000E-01 0.0 


0. 1000E-01 


0.0 


13 


0.1000E-01 0.9537E-05 0.1000E-01 


0. 0 


15 


0.1000E-01 0.0 


0. 1000E-01 


0. 9468E-02 


17 


0.1000E-01 0.0 


0. 1000E-01 


0.0 


19 


0.1000E-01 0.0 


0. 1000E-01 


0. 0 


21 


0.1000E-01 0.0 


0. 1000E-01 


0.0 


23 


0.1000E-01 0.1001E-04 0.1000E-01 


0. 0 


25 


0.1000E-01 0.9410E-02 0. 1000E+01 


0.0 


26 


0.1000E-01 0.0 


0. 1000E + 01 


0.0 


27 


0.1000E+01 0.0 


0. 1000E-01 


0. 1002E-01 
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TABLE XX 



COHPOSITE DRIVZSHAFT DESIGN SOHSARY: NORHAL 



COMPOSITE DRIVSSHAFT OOTPOT 



NUMBER OF PLYS = '4 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = 3.0 



DIMENSIONS 
PLY THICKNESS 

INSIDE 

1 0. 10000E-01 

2 0. 10212E+00 

3 0. 10212E+00 

4 0. 10000E-01 



PERCENT DIAMETER 
0.15 81 1E + 01 
4.46 0.17011E+01 

45.54 0.19054E+01 

45.54 0.21095E+01 

4.45 0.21295E+01 



THETA 

0.0 

0.3U968E+02 

-0.34968E+02 

0.90000E+02 



STIFFNESS 
AE = 0. 13972E+08 
El = 0.631 16E+07 
GJ = 0.56217E+07 



LOADS : 

L '' T 

r*0.31513E+05 



M F H.P. RPM 

0.0 0.500000E+02 0.15000E+03 0.30000E+03 



LOAD CONDITION 1 

PLY EPL S.F. EPT S.F. 

1 0.35785E-05 100.00 0.0 130.00 

2 0.25105S-02 3.41 -0.25069E-02 7.02 

3 -0. 27745E-02 3.09 0.27781E-02 1.70 

4 0.18741E-08 100.00 0.35766E-05 100.00 

CRITICAL SPEED = 3.11715E+04 

MAXIMUM DEFLECTION = 0.34815E-01 

WEIGHT = 0 .902 03S+01 

VOLUME = 0 . 161 08E<-03 



E^LT 

0.47679E-02 
0. 18288E-02 
0.20319E-02 
-0.59688E-02 



S.F. 

3.86 

10.06 

9.06 

3.08 
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"■'ili! ij! : 



r&BLE XXI 



COHPOSITE DRIVESHAFT COPES STANDARD OOTPOT: LOGNORMAL 



OPTIMIZATION RESULTS 
OBJECTIVE FUNCTION 

GLOBAL LOCATION 3 FUNCTION VALUE 0.10281E+02 



DESIGN VARIABLES 





D. V. 


G LOB AL 


LOW EH 




UPPER 


ID 


NO. 


VAR. NO. 


30U ND 


VALUE 


BOUND 


1 


1 


11 


0. 1C300E-0 1 


0. 19474E-01 


0 .50000E+00 


2 


2 


12 


o.ioaooE-oi 


0. 92916E-01 


0.50000E+00 


3 


2 


13 


0. 10000S-01 


0. 92916E-01 


0.50000E+00 


4 


3 


14 


0.1000 OS-0 1 


0. 10000S-01 


0.50000E+00 


5 


4 


32 


0.0 


0. 41392S+02 


0.90000E+02 


6 


4 


33 


0,0 


-0. 41 392E+02 


-0.90000E+02 


7 


5 


2 


0,5000 OE+0 0 


0. 20465S+01 


0.50000S+01 



DESIGN 


CONSTRAINTS 












DETERMI NISTIC 








GLOBAL 


LOWER 


MEAN 


UPPER 


ID 


VAR. NO. 


BOUND 


VALUE 


BOUND 


1 


51 


-0 .85700E-02 


0. 37986E-05 


0. 85700E-02 


3 


52 


-0.17600E-01 


0.0 


0.47100E-02 


5 


53 


-0 .18400 E-01 


0.34484E-02 


0. 18400E-01 


7 


54 


-0 ,85 700 E-02 


0. 18651E-02 


0.85700E-02 


9 


55 


-0 . 17600 E-01 


-0. 1 3613E-02 


0.47100E-02 


1 1 


56 


-0 . 18400 E-01 


0. 46804E-03 


0. 18400E-01 


13 


57 


-0.85700E-02 


-0, 201 32S-02 


0, 85700E-02 


15 


58 


-0.17600E-01 


0.20170E-02 


0.47 100E-02 


17 


59 


-0.18400E-01 


0.51418E-03 


0, 18400E-01 


19 


60 


-0 .85 700E-02 


0. 12861E-08 


0. 85700E-02 


21 


6 1 


-0 , 17600E-01 


0. 37974E-05 


0. 47100E-02 


23 


62 


-0.18400S-01 


-0.40960E-02 


0. 18400E-01 


25 


63 


0.31 513E+05 


0. 54975E+05 


0.11000S+16 


26 


65 


0 .30000E+03 


0. 12564E+04 


0. 11000E+16 


27 


6 6 


-0 .11000 E+16 


0. 29994S-01 


0. 50000E-01 
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r&BLE XXII 



COapOSITE DRI7ESHAFT COPES PROBABILISTIC OOTPOT: LOGNORHAL 



PROBABILISTIC 





MEAN 




STANDARD 


COEFFICIENT 


ID 


VALUE 


DEVIATION 


OF VARIATION 


1 


0. 37986E-05 


0.5215E-06 


0.5215E-01 


3 


0. 0 




0.0 


0.0 


5 


0. 34U84E-02 


0.5257E-03 


0. 1524E+00 


7 


0. 1865 IE-02 


0.2842E-03 


0. 1524E+00 


9 


-0. 1861 3E -02 


0.2838E-03 


0. 1525E+00 


11 


0. 46804E-03 


0.2815E-03 


0.6015E+00 


13 


-0. 20132E-02 


0.3069E-03 


0. 1524E+00 


15 


0. 20170E-02 


0.3073E-03 


0. 1524E+00 


17 


0. 5141 8E-03 


0.3062E-03 


0 .5956E+00 


19 


0. 12861E-08 


0.3216E-03 


0. 1000E+02 


21 


0. 37974E-05 


0.3215E-03 


0. 1000E+02 


23 


-0. 40960E-02 


0.6256E-03 


0. 1527E+00 


25 


0. 54975E + 05 


0.7813E+04 


0. 1421 E+00 


26 


0. 12564E+04 


0. 1012E+03 


0. 8058E-01 


27 


0, 29994E-01 


0.5135E-02 


0. 1712E+00 




PROBABILITY 


OF FAILORE 


PROBABILITY 


OF FAILURE 




AT LOWER 


BOUND 


AT UPPER BOUND 


ID 


ALLOWED 


CALCULATED 


ALLOWED 


CALCULATED 


1 


0. 1000E-01 


0.0 


0. 1000E-01 


0.0 


3 


0. 1000E-0 1 


0.0 


0. 1000E-01 


0.0 


5 


0. 1000E-01 


0 .0 


0. 1000E-01 


0.0 


7 


0. 1000E-01 


0.0 


0. 1000E-01 


0.0 


9 


0. 1000E-01 


0.0 


0. 1000E-01 


0.0 


1 1 


0. 1000E-01 


0 .0 


0. 1000E-01 


0.0 


13 


0. 1000E-01 


0 .0 


0. 1000E-01 


0.0 


15 


0. 1000E-01 


0.0 


0. 1000E-01 


0.3338E-04 


17 


0. 1000E-01 


0.0 


0. 1000E-01 


0.0 


19 


0. 1000E-01 


0.0 


0. 1000E-01 


0.5313E-02 


21 


0. 1000E-01 


0.0 


0. 1000E-01 


0.9957E-02 


23 


0. 1000B-01 


0.0 


0. 1000E-01 


0.0 


25 


0. 1000E-0 1 


0.1000E-01 


0. 1000E+01 


0.0 


26 


0.1000E-01 


0.0 


0. 1000E+01 


0.0 


27 


0. 1000E+01 


0 .0 


0. 1000E-01 


0. 1544E-02 
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TiBLE XXIII 



COHPOSITE DRIVESHAFI DESIGN SOHMAHY: LOGNORMAL 



COMPOSITE DRIVESKAFT OUTPOT 



NUMBER OF PLYS = ^ 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = 0.0 



DI MEN SIO NS 

PlV thickness percent DIAMETER 
INSIDE 0.20465E+01 
0. 19474E-01 9.04 0.20855E+01 

0.92916E-01 43.16 0.22713E+01 

0.92916E-01 43.16 0.24571E+01 



1 

2 

3 

4 



0. 10000E-01 



4.54 



0.24771E+01 



THETA 

0.0 

0.41392E+02 
-0.41392E+02 
0. 90000E+02 



STIFFNESS 
AE = 0. 13163E+08 
El = 0.82517E+07 
GJ = 0.95289E+07 



LOADS ; 

L ^ T 

r*0.31513E+05 



M F H.P. RPM 

0.0 0.5000C0E+02 0.15000E+03 0.30000E+03 



LOAD CONDITION 1 

PLY EPL S.P. EPT S.F. 

1 0.37986E-05 100.00 0.0 100.00 

2 0.18651E-02 4.59 -0.18613E-02 9.46 

3 -0. 20132E-02 4,26 0.20170E-02 2.34 

4 0.12861E-08 100,00 0.37974E-05 100.00 

CRITICAL SPEED = 0,12564B+0'4 

MAXIMUM DEFLECTION = 0.29994E-01 

WEIGHT = 0.10281E+02 

VOLUME = 0.18359E+03 



EPLT S.F. 

0.34484E-C2 5.34 

0.46804E-03 39.31 
0.51418E-03 35.79 
-0.40960E-02 4.49 
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T4BLE XXIV 



COMPOSITE DRIVESHiFT COPES OOTPOT: DETEEMI NISTIC 

OPTIMIZATION RESULTS 
OBJECTIVE FUNCTION 

GLOBAL LOCATION 3 FUNCTION VALUE 0.10231E+02 



DESIGN VARIABLES 





D. V. 


GLOBAL 


LOW ER 




UPPER 


ID 


NO. 


VAR, NO. 


BOUND 


VALUE 


BOUND 


1 


1 


11 


0.10000E-01 


0. 10000E-01 


0.50000E+00 


2 


2 


12 


0. 10000E-01 


0. 10697E+00 


0.50000E+00 


3 


2 


13 


0. 10000E-01 


0. 1 0697E+00 


0.50000E+00 


U 


3 


14 


0.1030 OE-01 


0. 10000 E-01 


0.50000E+00 


5 


4 


32 


0.0 


0. 30524E+02 


0 .90000E+02 


6 


4 


33 


0.0 


-0. 30524E+02 


-0.90000S+02 


7 


5 


2 


0.5000 OE+00 


0. 18376E+01 


0 .50000E+01 



DESIGN 


CONSTRAINTS 












DETERMI NISTIC 








GLOBAL 


LOWER 


MEAN 


UPPER 


ID 


VAR. NO. 


BOUND 


VALUE 


BOUND 


1 


51 


-0 ,42850 E-02 


0.26961S-05 


0. 42850E-02 


3 


52 


-0 .88000E-02 


0.0 


0.23550E-02 


5 


53 


-0.92000 E-02 


0.43686E-02 


0. 92000E-02 


7 


54 


-0 .42 850 E-02 


0.21334S-02 


0.42850E-02 


9 


55 


-0 .88000 E-02 


-0.21307E-02 


0. 23550E-02 


11 


56 


-0.92000E-02 


0.23560E-02 


0. 92000E-02 


13 


57 


-0.42 850 E-02 


-0.23496E-02 


0. 42850E-02 


15 


58 


-0 .88000 E-02 


0. 23522E-02 


0. 23550E-02 


17 


59 


-0.92000 E-02 


0.26043E-02 


0. 92000E-02 


19 


60 


-0 .42850E-02 


0. 17024E-08 


0. 42850E-02 


21 


6 1 


-0.88000 E-02 


0.26944E-05 


0. 23550E-02 


23 


62 


-0.92000E- 02 


-0.54219E-02 


0. 92000E-02 


25 


63 


0.63025 E+05 


0.54108E+05 


0. 11 OOOE+16 


26 


65 


0 .60000 E+03 


0. 1 3824E+04 


0. 11000E+16 


27 


66 


-0.11000 E+16 


0.24515E-01 


0. 25000E-01 
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riBLE XXV 



COHPOSITE DBIVESHAPT DESIGN SOMSARY: DETERMINISTIC 



COMPOSITE DRIVSSHAFT OOTPai 



NDHBER OF PLYS = % 

NOMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = 3.0 



DIMENSIONS 

PLY THICKNESS PERCENT DIAMETER 
INSIDE 0.18376E+01 

1 0. 10000E-01 a. 27 0.18576E+01 

2 0. 10697E+00 45.73 0.20715E+01 

3 0. 10697E + 00 45.73 0. 22855E+01 

4 0. 10000E-01 4.27 0.23055E+01 



THETA 

0.0 

0. 33524E + 02 
-0.30524E+02 
0. 90000E + 02 



STIFFNESS 
AE = 0.1 8546E+08 
El = 0.99273E+07 







GJ = 0.66999E+07 












LOADS : 










L.C 




T 


M 


F H.P. 


RPM 




1 


0 


.31513E+05 


o 

• 

o 


0.500000E+02 3.15000E+03 0.30000E+03 


LOAD 


CONDITION 


1 








PLY 




EPL 


S.F. 


EFT S.F. 


EPLT 


S.F. 


1 


0. 


26961E-05 


100. 00 


0.0 130.00 


0.43586E-02 


4.21 


2 


0. 


21334E-02 


4. 02 


-0.21307E-02 8.26 


0. 23560E-02 


7.81 


3 - 


0. 


23496 E-02 


3.65 


0.23522E-02 2.00 


0.26043E-02 


7.07 


4 


0. 


17024E-08 


100.00 


0.26944E-05 130.00 


-0.54219S-02 


3.39 






CRITICAL SPEED 


= 0 . 13824E<-04 










MAXIMUM DEFLECTION = 0.24515E-01 










WEIGHT 




= 0 .1023 1E+02 










VOLUME 




= 0.18270E+03 
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VII. CONCLUSIONS 



Numerical optimization using probabilistic design tech- 
niques provides an effective method for designing structures 
based on allowed component reliability. 

Existing ANALIZ subroutines can be easily modified to 
perform probabilistic optimization and as demonstrated, 
complex designs can be accomplished. 

The major drawback to the program developed is the large 
increase in function evaluations as compared to the standard 
deterministic optimization. This disadvantage can be 
partially offset by reducing the frequency in which the 
constrained variable gradients are calculated. 
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VIII. RECOMMENDATIONS 



A method of scaling the variables should be investigated 
to reduce ill-conditioning in the optimization process. 

This is a possible area to pursue in order to alleviate 
problems caused by very small coefficients of variation. 

Further efforts should be undertaken to accomplish the 
reduction of function evaluations while maintaining suitable 
optimization results. The altering of the frequency in which 
constrained variable gradients are calculated appears to be 
the most promising approach. 

In order to reduce function evaluations it is recommended 
that COPES be modified to allow the user to supply mathemati- 
cal expressions for the partial derivatives of the constrained 
variables via the ANALIZ subroutine. These could be used 
instead of the finite different gradients calculated by 
COPE19/COPE20 . Modifications to allow the user to supply 
precalculated gradients would also be helpful to this end. 

The sophisticated user can accomplish this now by using CONMIN 
directly in conjunction with C0PE19, C0PE21, and COPE22. 

It is further recommended that additional work be under- 
taken to modify COPES to allow for the use of the following 
options with probabilistic design: A) two variable function 

space; B) sensitivity analysis; C) optimum sensitivity. 

An additional worthwhile modification would be to provide 
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for the calculation in C0PE19 and subsequent output of the 
standard deviation of the objective function mean value. 

This investigation dealt only with component reliability 
in the design process. Work should be undertaken to ascer- 
tain the potential for the inclusion of system reliability 
considerations in the optimization process. 
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APPENDIX A 



COPES MANUAL ADDENDUM 

This appendix is intended as an addendum to the COPES 
Manual, [Ref. 1]. 



I . INTRODUCTION 

The purpose of this document is to provide user instruc- 
tions for performing numerical optimization using probabilis- 
tic design techniques. Four subroutines have been included 
in COPES (Control Program for Engineering Synthesis) to 
provide this additional design tool. 

This discussion describes the capabilities of probabilis- 
tic optimization using the COPES/CONMIN program. A simple 
design example is first presented to demonstrate the program 
capabilities. Guidelines are given for writing analysis 
codes. The data organization is outlined and sample data is 
presented. 

This publication was written to serve as an addendum to 
the COPES Manual; "COPES — A Fortran Program for Engineering 
Synthesis," L. E. Madsen and G. N. Vanderplaats , NPS69-81- 
003, Naval Postgraduate School, Monterey, California, March 
1982. 
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II. DESIGN EXAMPLE 



It is required to design the cantilevered beam shown in 
Figure 1. The objective is to find the minimum mean volume 
of material which will support the concentrated load and 
maintain an allowed probability of failure of 10%. 



P 




That is , 



Minimize volume = B * H * L (1) 

The mean bending stress in the beam must not exceed its 
limit which is 20,000 psi which has a coefficient of variation 
of 0.07; 



BSTRES 



Me 

I 




20,000 



( 2 ) 
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The mean shear stress must not exceed its limit of 10,000 
psi which has a coefficient of variation of 0.05; 



SHRSTR 



3P _ 3P 
2A 2BH - 



10,000 



(3) 



and the deflection under the load must not exceed one inch, 
coefficient of variation equals 0.02; 



DELTA 



PL 

3EI 



4PL^ ^ 

3 _ 

EBH 



1.0 



(4) 



Additionally, geometric limits are imposed on the mean 
dimensions so that; 



0.5 1 B < 5.0 (5) 

1.0 < H < 20.0 (6) 

I 1 10.0 (7) 

The manufacturing procedure is such that all dimensions have 
a coefficient of variation of 0.01. The mean value of the 
dead load is 10,000 lb. with a coefficient of variation of 
0.07. The Young's modulus mean value is 30.E+06 psi with a 
0.06 coefficient of variation. 

The ANALIZ subroutine on page 11 of the COPES manual 
is used for the analysis. Only one modification is needed. 
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The variable P, for the load, is added to the global common 
block. 

The COPES data used for the standard optimization, pages 
60 and 61 of the COPES manual is modified to perform proba- 
bilistic optimization. NCALC in data block B is changed to 
7 to indicate probabilistic optimization is to be performed 
and IPROB equals 1 is entered in column eight so the normal 
distribution model will be used. Three additional data 
blocks are required, Ul, U2 and U3. 

DATA BLOCK Ul : 

There are five probability variables which make up the 
constraint equations 2, 3, 4 and 7. They are B, H, E, AL, 
and P; so IVAR equals 5. It is desired that gradients of 
the constrained variables only be calculated at the beginning 
of each CONMIN iteration, therefore LGRAD equals 1 is entered 
in column two. 

$ DATA BLOCK Ul 
5,1 

DATA BLOCK U2 : 

The global location of the probability variables and 
their corresponding coefficients of variation are entered 
in this block. 

$ DATA BLOCK U2 



.01 


width, B 


1 — 1 
O 

• 


height , 
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10, .05 



load, P 

Young's mod, E 
length , AL 



8, .06 
9, .01 

DATA BLOCK U3 : 

This block contains the allowed probability of failure 
and coefficient of variation at both the upper and lower 
bound of each constrained variable. All four constrained 
variables have upper boxonds , but are unbounded on the lower 
end. 



$ DATA BLOCK U3 




1. ,.07, .1,.07 


BSTRES 


1.,.05,.1,.05 


SHRSTR 


1. , .02, .1, .02 


DELTA 


1 — 1 
o 

t 

» — 1 

• 

1 — 1 
o 

« 

• 

» — 1 


H/B 



The allowed probability of failure at the upper bound and 
the coefficient of variation are listed in columns three 
and four. Because there are no lower bounds the resulting 
probability of failure will be zero. These first two columns 
will be ignored during probability of failure calculations 
in COPES. 

Coupling the ANALIZ subroutine to COPES yields the 
following optimum design; 
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CANTILEVERED BEAM 



AL 


= 


200.0 


P 


= 


0 . lOOOOE+05 


E 


= 


0.30000E+08 


B 


= 


1.92 


H 


= 


18.75 


VOL 


= 


7261.54 


BSTRES 




0.17802E+05 


SHRSTR 


= 


0 . 41727E+03 


DELTA 


= 


0 . 84382E+00 


H/B 


= 


9.78 



The COPES output for this example is provided in Figure 2. 

Ill . PROGRAMMING GUIDELINES 

The programming guidelines presented in Chapter II of 
the COPES manual all apply in performing probabilistic 
optimization. One additional rule must be followed when 
writing the analysis siobroutine: All variables in the 

constraint equations which have probabilistic distributions 
must be included in the global common block of the analysis 
subroutine. 



IV. DATA BLOCKS 



Input instructions for the probabilistic optimization 
data blocks are presented in the following five pages. 
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of blocks A-K and V is required. 
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DESCRIPTION: Probability Variables and Gradient information. 
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bound coefficient of variation. Enter .0 if corresponding 
data block I is greater than l.OE+15. 



V. NOTES TO THE USER 



1. The probabilistic optimization feature is not 
compatible with the following COPES options; 

A. sensitivity analysis 

B. two variable function space 

C. optimum sensitivity 

D. approximate optimization 

2. It is not recommended that one attempt to design for 
allowed probabilities of failure of less than 0.00001 or 
greater than .99999 due to modifications made to "widen” the 
cxamulative normal distribution curve. 

3. The optimization process is very sensitive to the 
coefficients of variation of the probability variables which 
make up the constraint equations. Coefficients of variation 
of less than 0.0001 may result in erroneous results if the 
coefficients of variation of some of the constraint limits 
are also less than 0.0001. 

4. In some cases the initial design will be such that 
the optimization process terminates unsatisfactorily due to 
the calculation of zero constraint gradients. Termination 
will occur either due to an inability to achieve a feasible 
design or CONMIN will reach the maximum number of iterations 
no matter how large ITMAX is made. This probably will be 
due to the fact that the probabilities for failure calcu- 
lated are either 1.0 or 0.0 and the finite difference steps 
are insufficient to reach the region along the cumulative 
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normal distribution curve between Pf = 1.0 and Pf = 0.0. 

This could be caused by an initial design which is a great 
distance from the optimum or by very small coefficients of 
variation. If this happens, perform a deterministic optimi- 
zation to achieve a design in the vicinity in the design 
space desired and use this result for the initial design in 
probabilistic optimization. 

5. In order to determine probabilities of failure, COPES 
calculates the gradients of the constrained variables with 
respect to the probability variables which make up the 
constraint equations. This is done by finite difference 
steps resulting in numerous calls to the ANALIZ subroutine. 

In many problems these gradients remain essentially the same 
so the user is given the ability to determine the frequency 
in which these gradients are calculated. This provides for 
a significant reduction in the number of function evaluations. 
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Figure 
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nHleverec beam analysis and design 
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TOTAL NUMBER OF CONSTRAINED PARAMETERS 



INPUT CATA FGR PROBABILISTIC DESIGN 
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CANTILEVERED BEAM 
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FORTRAN PROGRAM FOR 
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